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ABSTRACT 

We present a study of the spatial and spectral evolution of the loop-top (LT) 
sources in a sample of 6 flares near the solar limb observed by RHESSI. A dis- 
tinct coronal source, which we identify as the LT source, was seen in each of these 
flares from the early "pre-heating" phase through the late decay phase. Spectral 
analyses reveal an evident steep power-law component in the pre-heating and 
impulsive phases, suggesting that the particle acceleration starts upon the onset 
of the flares. In the late decay phase the LT source has a thermal spectrum and 
appears to be confined within a small region near the top of the flare loop, and 
does not spread throughout the loop, as is observed at lower energies. The total 
energy of this source decreases usually faster than expected from the radiative 
cooling but much slower than that due to the classical Spitzer conductive cool- 
ing along the flare loop. These results indicate the presence of a distinct LT 
region, where the thermal conductivity is suppressed significantly and/or there 
is a continuous energy input. We suggest that plasma wave turbulence could 
play important roles in both heating the plasma and suppressing the conduction 
during the decay phase of solar flares. With a simple quasi-steady loop model we 
show that the energy input in the gradual phase can be comparable to that in 
the impulsive phase and demonstrate how the observed cooling and confinement 
of the LT source can be used to constrain the wave-particle interaction. 

Subject headings: conduction — turbulence — Sun: flares — Sun: X-rays — Sun: loop- 
top source 
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1. INTRODUCTION 

The properties of thermal and nonthermal emissions of solar flares and their correlation 
are very important in understanding the flare energizing process. The well-known Neupert 
effect (Neupert f 968) states that the time integral of the hard X-ray (HXR) or microwave 
emission is correlated with the light curve of the soft X-ray (SXR) emission. This effect 
is generally explained with the non-thermal thick-target model, where flares are triggered 
by injections of nonthermal electrons that produce the HXRs and deposit most of their 
energy at the footpoints (FPs). This causes chromospheric evaporation filling the loop with 
hot thermal plasmas responsible for the observed SXRs (Brown 1971; Lin & Hudson 1971; 
Petrosian 1973; McTiernan et al. 1999). However, detailed studies (Dennis & Zarro 1993; 
Veronig et al. 2002) show deviations from this simple effect, indicating that this picture is 
not complete. There is clear evidence of a SXR rise before the impulsive HXR emission of 
some flares (often referred to as the pre-heating) , and a persistence of the increase or slower 
than expected decrease of the SXR flux is frequently observed after the HXR emission has 
died out. Statistical studies (e.g., Lee, Petrosian & McTiernan 1995) indicate a heating 
process other than the coUisional heating mentioned above. 

Since the discovery of impulsive HXR emission from the loop-top (LT) by Yohkoh (Ma- 
suda et al. 1994; Masuda 1994), it has been realized that the spatial and spectral evolution 
of the LT source plays a crucial role in revealing the energy release process of magnetic recon- 
nection and the consequent particle acceleration and plasma heating in solar flares (Tsuneta 
1996; Alexander & Metcalf 1997; Petrosian, Donaghy & McTiernan 2002). To explain these 
observations one needs to trap high-energy electrons in the LT region. However, due to the 
low energy resolution of the Yohkoh HXR Telescope (HXT), the exact nature of the LT and 
FP sources is not well-constrained. It is, e.g., difficult to tell whether the LT source has 
a steep nonthermal spectrum or is produced by a "superhot" thermal gas (Lin & Schwartz 
1987; Tsuneta et al. 1997; Petrosian et al. 2002). Analyses of unresolved observations during 
the late decay phase have shown that the observed energy dechne does not agree with that 
expected through conduction (Moore et al 1980; McTiernan et al. 1993). Several mecha- 
nisms may alleviate this discrepancy (Antiochos & Sturrock 1976, 1978; Rosner et al. 1986; 
Takahashi & Watanabe 2000; Reeves & Warren 2002), however, these models were designed 
to explain observations with low spectral, spatial and temporal resolutions. RHESSI with 
its higher resolutions can be very helpful in this regard. 

In this paper we present a systematic imaging and spectroscopy investigation of the 
LT source of 6 simple RHESSI limb fiares, each of which seems to be associated with a 
single loop, instead of multiple loops for complex fiares. In general, the spectrum of the LT 
source has two components, a low-energy thermal component plus a power-law tail (or a 
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broken power-law). The nonthermal component dominates in the early phases. The thermal 
component becomes more and more prominent as the flare proceeds, and in the late decay 
phase the nonthermal component vanishes completely. From spectral fits we determine the 
evolution of the emission measure EM and temperature T of the LT source during the decay 
phase. We also determine the evolution of the volume of the source V . These give us the 
evolution of the density n and the total energy £ = SnkBTV, where /cb is the Boltzmann 
constant. Prom these we find that the energy decay rate is higher than the radiative cooling 
rate, but much lower than that due to the thermal conduction, consistent with previous 
studies (e.g., McTiernan et al. 1993). Clearly suppression of conduction and/or continuous 
energy release are required to explain these observations. Another important and surprising 
feature of these sources is that they are often found localized near the LT, which is distinct 
from observations at lower energies where emission along the whole length of the loop is seen 
in an arcade structure. We explore the possible mechanisms for the confinement of these 
sources and for the suppression of conductivity and show that the former alone rules out 
quasi-steady loop models with a uniform conductivity even when the radiative cooling is 
taken into account. We suggest that plasma waves or turbulence, which we believe plays a 
major role in the acceleration of particles during the impulsive phase (see e.g., Petrosian et 
al. 2002), can be the agent of both the heating and confining of the LT source during the 
decay phase as well. We consider a steady-state isobaric loop model with turbulence plasma 
waves concentrated in a distinct LT region and explore the degree to which the observed LT 
confinement can set constraints on the plasma heating and conduction suppression processes. 
The model predicts strong line emissions from the FPs in the gradual phase, and for flares 
with a relatively longer decay time the energy input in the decay phase may exceed that of 
the impulsive phase. This is in direct contradiction with the Neupert effect, which claims 
that the energy release is limited to the impulsive phase. Broadband observations, especially 
those in the optical and UV bands, should be able to test the model. 

In § 2 we present the analyses of these flares and the key observational features. The- 
oretical investigations of these results and their implications are discussed in § 3, where a 
viable model is also presented. § 4 summarizes the main conclusions of this study. 

2. OBSERVATION AND DATA ANALYSIS 

In this section we flrst discuss the flare selection criteria. The spatial and spectral 
properties of six selected flares are then studied in detail. 
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2.1. Flcire Selection 

Because we are interested in the evolution of flares, a good time-coverage of the whole 
event is required. To reduce the instrumental effects, only events with their shutter state 
remaining unchanged are considered. We also focus on relatively large flares with peak count 
rate > 300 per detector in the 6 — 12 keV energy band to ensure a high image quality 
and rehable spectral results. To study the spatial evolution of the LT source, the flares need 
to be close to the solar limb and have a relatively simple morphology. Furthermore, to avoid 
spectral contaminations from the FP sources, flares with their LT source staying above the 
solar limb and their FPs completely occulted by the solar disk would be ideal candidates. 
We choose flares with heliocentric longitude larger than 75°. 

We searched events observed by RHESSI in 2002 and found a total of six flares appro- 
priate for this study. Two of them are limb flares with both their LT and two FP sources 
clearly seen, one is a partially occulted flare, where the LT source (above the limb) and one 
FP source were observed, and the three remaining flares are coronal events above the solar 
limb (see Figures 1-3). For each flare we also examined the simultaneous G'Oi?^ observations 
extending to lower energies, and SOHO EUV Imaging Telescope (EIT) 195 A observations 
to determine the loop geometry. From these and observations described below, we found 
that all of these flares appear to be associated with a single-loop structure. 

2.2. Images 

We obtained images of the flares with the PIXON algorithm, which has the best spatial 
resolution. Because the energy resolution of detector 2 is severely degraded, which introduces 
signiflcant uncertainties to the source structure at a given energy range, only the front 
segment of detectors 3-8 were used, achieving a resolution of ~ 5'/0 (Aschwanden et al. 
2003). These images reveal a distinct LT source present from the beginning of the flares 
through the late decay phase, whereas the FP sources, whenever observed, are seen only 
during the impulsive phase at higher energy HXRs. Figures 1-3 show the 6-12 keV images 
of the six flares at the HXR peak (thin contours for the first three and grey scale for the 
rest) and a late period in the decay phase (thick contours). The peak time FP sources (at 
51-57 keV, 34-39 keV, and 60-100 keV, respectively) are indicated by the gray scale for the 
two hmb flares and the partially occulted April 4 flare^. The plus signs indicate the emission 



^Note that one of the high energy sources for this flare, which would be normally associated with the 
FPs, is above the solar limb. 
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centroid of the LTs. There are two flares on April 4, one of which is partially occulted labeled 
"a" and the other is a corona event labeled "b" (Figure 2). For the April 30 flare, there 
are two corona sources during the impulsive phase. The lower one is associated with the 
LT, while the higher one may indicate the other end of a rcconnection current sheet (Sui & 
Holman 2003). The dotted contour at the upper-right corner gives the FWHM of the PSF 
of each flares at 6-12 keV. We see that all of the LT sources are partially resolved. 

Although the LT sources appear to shift systematically upward as the flares proceed, the 
motions are not significant. And the LT sources seem to be confined within a small region 
with no significant changes in the source morphology throughout the fiarc. These suggest 
that the flares may be associated with a single-loop. More importantly, the confinement 
of the image indicates that they are not in line with the source evolution expected in the 
thick-target chromospheric evaporation model (Mariska et al. 1989). On the other hand, 
all the flares studied here underwent a prominent preheating phase. It is possible that this 
preheating process raised the gas pressure in the flare loop and suppressed the post impulsive 
phase chromospheric evaporation effects (Emslie et al. 1992). The results are also distinct 
from flare images obtained at lower energies, such as those taken with the KIT and Yohkoh 
SXT (1 — 2 keV), where the complete arcade loop structure is often seen. 

For quantitative analyses we assume that the LT source is an ellipsoid and determine 
its semi-minor (perpendicular to the loop) and semi-major (along the loop) axes a and b by 
fltting the 15% contour of the observed peak emission (outer contours in Figures 1-3 and 
red outer contours in Figure 4). To compare the observed source image with the theoretical 
loop models discussed below, we assume that the loop is a half circle and slightly tilted with 
respect to the solar surface so that the apex of the loop falls at the observed LT position, 
and the two FPs are at the observed FP emission centroids. The black fines in Figure 4 
show this loop structure (curved) and the direction perpendicular to it (straight). In Figure 
5 the brightness profiles along (upper) and in perpendicular to (lower) the loop are shown 
by the solid lines, the dot-dashed lines give the PSF. We see that the September 20 fiare is 
resolved. The August 12 fiare is resolved along the loop, but not resolved in the perpendicular 
direction. Because of this the size b and specially a should be considered as upper limits so 
that the volumes (densities) deduced below are upper (lower) limits. 



2.3. Total Spectrum and Imaging Spectroscopy 

For flares with their FPs completely occulted by the solar disk, the spectral evolution of 
the LT source over the entire flare duration can be readily studied with the spectral software 
package SPEX. The dominance of the LT source during the rising and decay phases also 
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facilitates studying the spectrum of the LT source without employing the comphcated imag- 
ing spectroscopy algorithms for partially occulted or unoccultcd flares. Because quantitative 
imaging spectroscopy of solar flares are limited by the low dynamical range, we only consider 
the imaging spectroscopy at the HXR peak of the two limb flares with both FPs seen and 
the evolution of the total spectrum, assuming that the impulsive phase spectral evolution of 
the LT source of the three flares with FP sources is similar to that of the LT source of the 
three occulted flares. 

For the August 12 flare we studied the spectra above 8 keV to avoid emission line 
features at ~ 7 keV. SPEX can handle these emissions in a thermal model but not for the 
nonthermal ones. For the rest of the flares we limited the spectral fltting to the above 10 
keV energy range because at least one of the shutters was in for these flares, which results in 
low count rates and a non-diagonal detector response matrix at low energies (below 12 keV). 
This causes significant uncertainties in the spectral fitting. In many cases, both a thermal 
and a power-law model can fit the observations with similar reduced . 

The left panel of Figures 6-8 show the best-fit spectra of the two hmb fiares of September 
20 and August 12, and the occulted fiare b of April 4 at the preheating phase (lower) and 
HXR peak (upper). For the hmb fiares, the high-energy HXR emission in the impulsive phase 
mostly comes from the FPs. The right panels are for an early (upper) and a late (lower) 
period in the gradual decay phase. The bottom panel of each figure gives the residual for the 
upper spectrum in unit of the standard deviation. These four periods arc indicated by the 
arrows in the top panels of Figures 10-11. With RHESSFs superior spectral resolution and 
broad energy coverage, it can be shown clearly that although the power-law component can 
still be seen early in the decay phase, the spectrum becomes purely thermal later on. The 
observations then reveal a hot plasma confined at the LT, and the confinement is evident at 
the highest energy. It is obvious that a combined study of images at different energy ranges 
will disclose the temperature and density profiles of the fiare loop in more detail, providing 
critical constraints on its dynamical evolution, especially the heating and cooling processes. 
Such kind of investigations are clearly warranted. 

Figure 9 shows the LT and FP spectra of the September 20 (left) and August 12 (right) 
fiares at their HXR peak. For the former only the thermal component of the LT can be seen, 
while the thermal component of the FP spectra are likely due to contaminations by the LT 
considering the compactness of the source and the limitation of the spatial resolution of the 
PIXON algorithm. Fitting the high energy spectra of the FPs with a power-law model, we 
obtained a spectral index of 2.75 and 2.65 for the southern and the northern FP, respectively, 
indicating an injection of a power-law electron spectrum into the FPs. The bremsstrahlung 
yield of this thick target source is then well defined (see e.g. Petrosian 1973), and gives an 
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electron energy flux of ~ 2.4 x 10'^'^{EjlOkeV)-^-^^ and ~ 0.8 x 10^'^{EjlOkeV)-^-^^eTg s'^ 
for the two FPs, where is the low energy cutoff of the power-law electron distribution. 
These are much higher than the bolometric luminosity of the thermal LT source: ~ lO^^erg 
s~^. For the August 12 flare the LT spectrum can be determined only below 10 kcV and 
therefore is not shown here. The electron fluxes at the southern and northern FPs are ~ 
2.4 X 10^^{Ec/WkeV)-^-^ and ~ 2.4 x 10^^{Ec/lOkeV)-^-'^erg s-^ respectively. The impulsive 
phase durations of the September 20 and August 12 flares are ~ 100 and ~ 30 s, giving a 
total energy input of ~ 3.2 x 10^^ and ~ 1.4 x 10^^ erg, respectively. 

2.4. Spectral Evolution 

The evolution of a flare is usually divided into three phases: "the prc-hcating phase", 
when the SXR emission rises monotonically without apparent HXR emission; the impulsive 
phase, when almost all of the high-energy (> 25 keV) HXR emission is produced, and the 
emission shows short time scale (< 1 sec) variations; and the decay phase, when the HXR 
emission vanishes and the SXR emission decays gradually. For some flares the pre-heating 
phase may not be obvious or be totally absent. The impulsive and decay phases, however, 
are common features of all flares. The top panels in Figures 10-12 give the light curves of 
the six flares at 6-12 keV, 25-50 keV and 50-100 keV, showing these three phases. 

As expected the spectral properties of a flare during these three phases are quite differ- 
ent. To study the spectral evolution, wc flt the total spectrum of these flares with a thermal 
plus a high-energy power-law or a broken power-law tail (the power-law index is usually 
larger at higher energies). The evolution of the power-law spectral index(es) 6 (dot-dashed), 
emission measure EM (dashed), and temperature T (solid) are shown in the second panel of 
Figures 10-12. To determine the spectral parameters with a reasonable certainty, one needs 
at least 4 spectral data points for each of the thermal and nonthermal components. Only 
parameters, which are well measured, are shown in the flgures.^ 

In the preheating and early impulsive phases, there is no clear evidence for a thermal 
component above 10 keV; only the power-law component is seen. The power-law component 



^Because the power-law component for the LT source is always very steep (with a photon index larger than 
5), it will dominate the thermal component at very low energies during the period when both components 
are present. We consider this scenario physically unreasonable, because the nonthermal particles arc likely 
accelerated from the thermal background. The power-law component should be truncated above the thermal 
energy of the background plasma. We therefore introduce a low-energy turnover of ^ 16 keV for the power- 
law component (below which the photon spectrum is assumed to be flat) so that the thermal component 
dominates the low-energy emission. 
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vanishes early in the decay phase, and the thermal component often appears after the HXR 
peak and gradually becomes more and more dominant. The fact that the nonthermal com- 
ponent is present from the rising phase to the early decay phase suggests that the particle 
acceleration or energy release occurs over a period longer than the duration of the impul- 
sive phase. There also appears to be a correlation between the thermal and the power-law 
components. The emission measure, EM, of the thermal component keeps increasing as far 
as the nonthermal component exists in qualitative agreement with the thick target model 
and the Neupert effect. The temperature, T, however, starts to decay minutes before the 
nonthermal component vanishes, qualitatively consistent with predictions of the evaporative 
coohng model (Antiochos & Sturrock 1978). However, it is not clear whether this evapora- 
tion is driven predominately by the nonthermal component or by a heat flux associated with 
the thermal particles. More detailed analyses and modeling are required to address these 
quantitatively. 

In the three flares with detected FP sources there are no signs of FP emission minutes 
after the HXR peak, when a relatively steep nonthermal component is still present in the total 
spectrum, indicating continuous particle acceleration conflned to the LT region. Evidently, 
the electrons accelerated during this phase are stopped in the coronal loop by the dense 
thermal plasma formed after the HXR peak (Veronig & Brown 2004). 



The most interesting flnding of this study occurs in the late decay phase, when the flare 
spectra are fitted well with a simple thermal model, while images show that the LT source 
is confined within a small region and remains stable there for several minutes. 



From the images and spectral analyses we obtain the evolution of the volume V — 
47ra^6/3, EM and T of the LT source. Here following the standard procedure (see McTiernan 
et al. 1993) we have assumed that the line-of-sight semi-axis of the source is equal to a. For 
a volume filling factor equal to unity (see, however, below), one can get the mean electron 
density rig, the total energy and the pressure P of the thermal LT source with the following 
equations: 



3. 



THEORETICAL MODELING 



3.1. Cooling Processes: Theory vs Observation 



£ = 



{EM/Vf^ = {'iEM/A-Ka%Y/'^ , 
?,nekBTV = ?,EM kBT/ue , 



(1) 
(2) 
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P = 2nekBT, (3) 

where, for the sake of simphcity, we have assumed a fully ionized pure hydrogen plasma 
(so that the proton density Up = Ug = n) and that the electron and proton temperatures 
are equal. One then obtains the density and energy evolution and the energy decay rate 
T^^ = where £{< 0) is the time derivative of E. The density evolution is indicated by 

the dotted hue in the second panel of Figures 10-12, and the energy decay rate by the crosses 
in the third panels. Note that for the flare of March 28, the energy increases occasionally in 
the gradual phase. The corresponding increase rate is indicated by the star signs. 

There are three major cooling processes: expansion, radiation, and conduction. As 
mentioned above there is no obvious expansion at the edge of the LT source, and adiabatic 
expansion down the legs of the arcade loop may not be important here since the pressure 
equilibrium is established on a time scale of the sound travel time r^, ~ L/cg ~ 20 s, where 
L ~ 10^ cm is the half-length of the loop and Cg ~ 400 km s^^ is the sound speed of the 
thermal plasma in the loop. This time scale is much shorter than the energy decay time 
~ 10 mins. 

For plasmas with a temperature below 3 keV, the radiative cooling is dominated by 
line emissions (Tandberg-Hanssen & Emslie 1988) ^, and the radiative cooling rate can be 
approximated as 

£rad = KrulT-^/^ with Kr = 1.42 X IQ-^'^erg cm^s-^K^/2 , (4) 

where all quantities (here and in what follows) are expressed in the cgs units unless specified 
otherwise. The corresponding energy decay rate 

-,-i,^V ,,„,„_4-i^ EM y- f T Y"'' /a ^-1 / b 



,,-.,^.,7xlO-V'^j55^J (5) 

is represented by the dash-dotted line in the third panel of Figures 10-12. (At the distance 
to the Sun, 1" ~ 7.3 x 10^ cm.) The radiative cooling rate is lower than the observed energy 
(or temperature) decay rate by a factor of a few, and the ratio of the two varies with time 
and is different for different flares, suggesting that the radiative cooling could not be the 
dominant cooling mechanism. This conclusion is different from that drawn by McTiernan et 
al. (1993), who used the instrument SXT that operates at lower energies and is most sensitive 
to relatively colder plasmas. The RHESSI LT source could be different from the SXT source 
because the former is usually hotter than the latter by ~ 10^ K. The fact that the RHESSI 
LT source is much smaller than that of the Yohkoh further supports this interpretation. 



^Clearly here we are not dealing with a pure hydrogen plasma. 



-10- 



The standard thermal conductive cooling rate is given by Spitzer (1962). For the LT 
sources understudy the conduction heat flux density along the loop and the corresponding 
cooling rate are 

^spit = KsT^^^VT with = 1.0 X lO^Vg cm-V^K-^/^ (6) 

-1 _ -^SpitTTO^ 



'^Spit 



where we have used the approximation VT" ~ T/L in obtaining the last expression. The 
thick solid hne in the third panel of Figures 10-12 show the energy decay rates expected from 
the Spitzer conductivity. As evident the conductive cooling rate is much higher than the 
observed energy decay rate. This result qualitatively agrees with the work by McTiernan et 
al. (1993). 

However, laboratory experiments and numerical studies have shown that the actual 
conductivity is smaller than that given by Spitzer when the temperature variation length 
scale L ~ 10^ cm is less than 30 times longer than the mean free path of the thermal electrons 
A = 1.4 X 10\T/WKy{n/10^^cm-^)-^ cm (Luciani, Mora, & Virmont 1983; Post 1956). 
This is the case for typical solar flare conditions. As a result the electrons deviate signiflcantly 
from a Maxwellian distribution, and the conduction is suppressed. This is referred to as the 
non-local conduction since the energy transport depends upon the global structure of the 
system. In principle, one needs to solve the Fokker- Planck equation to address this problem. 
Fortunately, the cooling rate can be approximated fairly well with the expression (Rosner, 
Low, & Holzer 1986) 

/ X N -0.36 

T„r^~ 0.11 Uj rspit-^ (8) 

The thin solid lines in the third panel of Figures 10-12 show the expected cooling rate based 
on equation (8), which is lower than the Spitzer rate but still higher than the observed energy 
decay rate. 

Takahashi & Watanabe (2000) observed similar discrepancy between the decay rate of 
Yohkoh flares and that expected from the Spitzer conductivity and introduced a volume 
fllling factor 1] < 1 to resolve it. For a loop consisting of many thin fllaments the combined 
cross section area a'^ = rjo^. Since r^j^. oc a^/(£'MV)^/^ oc a, this means Tgpj^. oc r]^^"^ 
and can be reduced for 77 <IC 1 or densities much greater than those calculated above. In 
this case, however, one must also consider the change in the radiative cooling rate oc 
{EM/Vy^'^ oc oc 77"^/^, which increases with a smaller filling factor. Thus the total 
expected coohng rate r^l — T^rj'^^^ +Tgp^^r]^^^ , which now has a minimum of 2(TspitTrad)~^''^ 
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at ricT — ''"s^it/'^rad- With more accurate temperature and emission measure obtained with 
RHESSI we can calculate this minimum for each flares. The ratio of these rates to the 
observed cooling rate t^^ and the corresponding rjcr are shown in Table 1 for a typical time 
in the gradual phase of the six flares. As evident the minimum total cooling rate is in 
general at least one order of magnitude higher than the observed rate, and a low (< 0.01) 
filling factor is required, implying a density on the order of 10^^ cm~^ and a gas pressure 
comparable to or even higher than the magnetic field pressure. Both of these are unlikely to 
be the case in solar flares. We therefore conclude that a continuous heating (well after the 
nonthermal component vanishes) and/or a suppression of the conduction in the LT source 
are required to explain the observed low energy decay rate. In what follow we set 77 = 1 
unless specified otherwise. 



3.2. Models with the Spitzer Conductivity 

As shown in the previous section, the flare decay time is much longer than the micro- 
scopic and hydrodynamic times, the plasma in the loop reaches an isobaric quasi-steady 
state. Without the suppression of conduction the LT source must be continuously heated at 
a rate comparable to the conduction rate to explain the observed energy decay. In this case, 
however, we expect a nearly uniform emission flux along the corona loop, which disagrees 
with the observations discussed in § 2.2, where the thermal source was found localized within 
a small region at the top of the flare loop. 

To quantify this study, we consider a loop model with a heat flux injection confined 
to the very top of the loop (not along the loop). The LT temperature T^t is determined 
from the spectral fitting, and those at the FPs are chosen to match the chromospheric value 
Tpp ~ 0.2 MK. (We do not attempt to include the even colder photosphere region because 
other effects, such as convective turbulence, will make our model inapplicable. See e.g. Spicer 
1979.) The magnitude of the heat flux clearly depends on the conduction. The gravitational 
effects can be neghgible because of the high gas pressure. With the radiative and conductive 
cooling processes, and the observed energy decay rate S taken into account, we have the 
energy conservation equation within the loop 

^^^^i^--^^M-S/V, (9) 

where I is the distance along the loop from the middle of the LT source, and JF^onfi gives 
the heat flux density due to conduction. Since the loop is in pressure equilibrium, £ is 
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independent of I. For loops with constant cross section (i.e., a independent of I), we get 



where S is the conduction suppression factor. Solving this equation with S = 1, wc obtain the 
temperature profile along the loop shown by the dashed lines in the upper panels of Figure 13 
for the flares on September 20 (Left panel: ne = 8 x lO^^cm'^, Tlt = 2.3 x 10^ K, L = 11.5") 
and August 12 (Right panel: rie = 4 x lO^^cm'^, Tlt = 1-9 x 10^ K, L = 19.5"), where the 
physical conditions are for the periods right after the power-law component vanishes. Since 
a nearly isothermal heating region will make the LT source even larger than the observed 
source size, we ignore the size of the heating region (see discussion below). As evident, in 
this case the expected radiation along the loop will be nearly uniform and not as highly 
confined to the LT region as observed in the RHESSI band. 

We are mostly interested in the brightness profile along the loop, and the above analyses 
give only this variation. The source structure in perpendicular to the loop (and the field line) 
is not constrained by these analyses. In order to compare the model with observations we 
convolve the model predicted emission from the one dimension loop with a two dimensional 
Gaussian profile with its width determined by the observed semi-minor axis a. We simulate 
such a loop structure as it would appear in the 6 — 12 keV band when observed with RHESSI 
and analyzed with the PIXON algorithm. The blue curves in Figure 4 show the 15% and 75% 
contours of the simulated images. It is evident that the model predicted emission extends 
further towards the FPs than the observed emission (red contours). Figure 5 compares the 
model predicted brightness profiles (dashed lines) with observations. The agreement in the 
perpendicular direction is a result of our convolution (lower panels). But along the loop 
(upper panels), which is what matters here, the LT source is extended with an FWHM 
~ 30% and 50% longer than that observed for the September 20 (left) and August 12 (right) 
flares, respectively. These arc longer than the observed source lengths by about one FWHM 
and one and a half FWHMs of the PSF (shown as dot-dashed profllcs), respectively. 

The lower panels of Figure 13 show the energy fluxes Fin and Font injected at the LT 
(plus) and the FPs (triangle), respectively, the total radiative loss |£^rad| = / £'raddl^ (cross), 
and \S\ (diamond). The energy balance is achieved here by introducing a high heating 
rate at the LT, most of which is conducted to the FPs. F-^n — -^out ~ 10^^, and 10^^ erg 

for the September 20 and August 12 fiares, respectively. These are almost two order 
of magnitude higher than the respective radiative loss rates. There are several difficulties 
with this scenario. The most important one is that this high rate (over several minutes) 
implies a total energy input > 10^° erg for the September 20 flare (with a duration > 100 
s), which is much larger than the energy input of ~ 3 x 10^^ erg estimated in § 2.3 for the 




(10) 
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impulsive phase. If this energy is converted into radiation at the FPs, there will be significant 
radiative signatures (Spicer 1979), which is not observed. Most of this energy therefore has 
to be deposited into the cold chromospheric plasma and would cause evaporation, which 
appears to contradict the observed slow decrease of the EM during this phase. One possible 
(but unlikely) scenario is that most of the plasma evaporated from the FPs has a very low 
temperature and does not contributed to the observed X-ray emission. Detailed analyses of 
GOES observations of these flares may shed hght on this problem (AUred et al. 2005). 

With a volume filling factor ^ 1, the conductive heat flux can be reduced. However, 
this will not alter the temperature proflle signiflcantly, so that the discrepancy with the 
conflnement of the LT remains unresolved. An upward flow of evaporated plasma from the 
FPs can reduce the conductivity, but it also will not modify the temperature profile of the 
loop significantly (Antiochos & Sturrock 1978). The recently proposed multiple-loop models 
for solar fiares may explain the discrepancies in the cooling rates, since the observed energy 
decay is controlled by the assumed loop formation rate (Reeves & Warren 2002; Warren & 
Docshek 2005). However, this model also can not explain the confinement of the LT source. 
Therefore all these effects may account for the discrepancy in the cooling rates and reduce 
the apparent size of the LT but fail to explain the observed strong conflnement due to the 
nearly uniform temperature proflle resulting from the sensitive dependence of the thermal 
conduction on the temperature. We conclude that in the late decay phase these observations 
suggest that the conductivity has to be suppressed signiflcantly. 



3.3. More Realistic Models with Suppressed Conductivity and Heating 

The above discrepancies between the observations and the simple model, primarily due 
to the conflnement of the LT source, can be overcome with a suppression of conductivity. 
As we will see below even in this case the LT source must be heated continuously (but at 
a slower rate). We propose that plasma wave and turbulence produced at the LT can be 
responsible for both these effects. In this section we investigate the characteristics of such a 
scenario. 

We flrst point out that a uniform suppression of conductivity along the whole loop 
cannot resolve these discrepancies. We have solved equation (10) for various values of the 
suppression factor 5" > 1 and found that the temperature proflle changes httle with S. This 
is because the conductive heating flux and radiative cooling rate are very sensitive to the 
temperature, the temperature has to be nearly uniform to satisfy the boundary conditions. 
When the thermal conduction dominates the loop must be nearly uniform to carry an almost 
constant heat flux for given temperatures at the LT and FPs. When the radiative cooling 
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becomes important, (erad t?T~^I'^ oc T~^/^ for the isobaric model and for the range of 
temperature of interest here), the coohng rate is also very sensitive to temperature, and the 
plasma is in an unstable state. The loop temperature again has to be nearly uniform to 
satisfy the boundary conditions at the FPs"^. 

Therefore we consider a model, where the suppression and heating (and the turbulence 
responsible for them) are confined to the LT region. Roughly speaking the loop consists of 
two parts: a hot LT region where the observed source locates and cold legs of comparable 
size. 

We have compared with observations many simulated images with different temperature 
profiles as an input and found that the temperature of the plasma near the 15% contour of the 
observed images has to be at least two times lower than the measured source temperature. 
We first note that as a consequence of the isobaric assumption this means that independent 
of the details of the model heating may be required because the radiative cooling rate will 
be at least 2^/^ ~ 6 times higher in the legs than in the hot LT source. For the September 20 
flare, this will exceed the observed energy decay rate. The conductive cooling rate (oc T^/^) 
will be > 6 times lower (here we assume that the mean temperature gradient for the loop 
does not change), but can still be the dominant cooling process, implying a suppression of 
conduction in the LT region by at least a factor of 6. Note that the ratio of conductive to 
radiative cooling rates is proportional to so that if the temperature of the relatively cold 
legs Ti^eg < TiT(Tspit/Trad)^^^ thc radiative cooling in the leg will exceed the available heating 
carried by conduction from the LT region. The requirement of T^g^ < 2Tlt then implies that 
Trad/Tspit Hiust be larger than ~ 32 for the loop model to be applicable. For the six flares 
understudy this requirement is satisfled. (It is marginal true for the partially occulted flare 
"a" on April 4). The weak dependence of this hmit on the ratio of the observed conductive 
and radiative cooling rates also implies that the temperature in the legs can not be more 
than ~ 3 times lower than the measured temperature of the LT source (Table 1). 

In summary we require a suppression of conductivity and heating in the LT source but 
the heating rate must be lower than the value obtained in the previous section, where the 
conduction is not suppressed. 

A continuous but slow production of plasma waves or turbulence via a relatively slow 



^In principle, one can prescribe a heating function along the loop so that the radiative cooling process 
can be balanced by the divergence of the conductive heat flux, the observed energy decay and the local 

heating (Rosncr, Tucker, & Vaiana 1978). A steep temperature profile in consistence with the observed LT 
confinement might be realized if the heating function were fine-tuned. We consider this scenario unphysical 
and will ignore it in the following discussion. 



- 15 - 



magnetic reconnection above the LT in the decay phase may explain the locahzed suppression 
of conduction and heating in the LT region. In such a scenario the conductivity is suppressed 
due to the scattering of high-energy particles in the hot plasma by the waves, which reduces 
the particle mean-free-path or scattering time and suppresses conduction. The dissipation 
of the wave energy, on the other hand, can result in heating of the plasma. Here we consider 
a simple case, where the scattering rate of particles by the waves is independent of the 
particle velocity v (and therefore the temperature T) but is proportional to the turbulence 
energy density. The corresponding conductive heating flux is given by (Spicer 1979): 

=ond — -7—. 1 — -^Spit = fi I ^cn\rr5/2^ ' 

1 + Tcoui/rsc \l + S {l)T^'^ J 

where TqI^^^ ~ 15(T/10^K)~^/^(n/10^°cm~^) s~^ is the mean Coulomb collision rate of the 
thermal electrons carrying the heat flux, and the ratio of the mean wave scattering to coulomb 
collision rates is given by t~^/tcou1 ~ 5'(/)T^/^ (assuming isobaric condition). The sup- 
pression factor S will be proportional to the wave energy density and will be expected to 
decrease toward the FPs. In what follows we describe the spatial variation of the sup- 
pression by a Gaussian function of width w centered at the centroid of the LT source, i.e. 
S — ySoCxp {—f'/xiP) T^^^"^ . If the waves and particles are coupled via a resonance process, 
the corresponding particle acceleration rate r'^ ~ ^{vp^jvYT'^ ^ where vp^ — B / •^jAmmip is 
the Alfven velocity {B is the magnetic fleld, rup proton mass), and the coefficient ^ depends 
on the wave spectrum and details of the wave-particle coupling (Schlickeiser 1989). The 
energy conservation equation (9) is thus modified to: 

d(a^.Fcond) . ^ , ^ < 7--^ > 

a^Z = -^-d(0 -y + ^ , (12) 

where < >, the average of the acceleration rate over the particle distribution, denotes 
the heating rate. For a loop with a constant cross section, this equation becomes 



d/ Vl + '5(/)T5/2 d/ 

where Vth = (3/cB7'/me)^^^ is the thermal velocity of the electrons and rrie is the electron 
mass. Because heating is distributed throughout the LT region, one can set the heat flux 
at the top of the loop, i.e. at / = 0, equal to zero, which gives a boundary condition. We 
keep the other boundary conditions of the LT and FP temperatures the same as those in 
the previous section. With such a model the temperature profile can be obtained once one 
specifies 5*0 and w. However, to satisfy the boundary conditions at the LT and FPs, the 
parameter ^ has to attain certain specific values, implying that heating is required. This is 
essentially an eigenvalue problem with ^ being the eigenvalue and T{1) the eigenf unction. 
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For the physical conditions used before for the flares on September 20 and August 
12, in Figure 13 the solid lines in the upper panels show several of these profiles with the 
corresponding 5*0 and w indicated in the caption. The lower panels show the dependence 
of the loop energetics on 5*0, where F^^ (thick solid), F^ut (thin solid), f^j-ad (dotted), and S 
(dashed) give the spatially integrated heating by the waves, the heating flux injected into the 
FPs, the spatially integrated radiative cooling, and the observed energy decay, respectively 
The dot-dashed line shows the dependence of ^B'^ on 5*0 with the scale indicated on the right 
frame. As expected Fin and Fout decrease with 5*0. These are consequences of the conduction 
suppression, which reduces not only the level of heating needed to keep the energy balance 
but also the conductive heat flux injected into the FPs. The parameter ^5^, therefore, 
decreases with 5*0. The radiative cooling rate, on the other hand, increases with 5*0 because 
the overall temperature of the loop decreases (see the top panel). As expected the observed 
energy decay rate only weakly depends on 5*0. 

It is clear that the energy balance is achieved by the interplays of the conduction, 
radiation, heating, and the observed energy decay. For the September 20 flare the heating 
process always dominates and is balanced by conduction and radiative cooling with the 
latter being more important with the increase of 5*0. The observed energy decay is slow 
and therefore not important. These are quite different from those of the August 12 flare, 
where the observed energy decay dominates at large values of -S'o, and for small values of 
Sq the conductive cooling dominates. Interestingly, for both flares Sq has a maximum, 
where Font — 0, indicating no chromospheric evaporation in the phase, which is in line 
with the slow changes or even decay of the observed EM (Fig. 10). Above this maximum, 
no physical solution exist because conduction is not efficient enough to balance the rapid 
radiative cooling near the FPs. The models with the maximum conduction suppression at 
the LT are indicated by the thick solid lines in the upper panels. The simulated images for 
these models are indicated by the green contours and dotted lines in Figures 4 and 5, which 
agree with the observed LT confinements. 

These models therefore account for the observed flare decay. Because there are no 
heat injections at the FPs, the radiative cooling process plays a dominant role then. To 
further illustrate the details of the model. Figure 14 shows the flare energetics along the 
loop. It is important to note that more than half of the radiation is produced near the FPs. 
Strong optical and UV line emissions are thus expected from the FPs. The model predicted 
bolometric luminosities (at the periods of interest) of the September 20 and August 12 flares 
are ~ lO^^erg and ~ 2 x lO^^erg s~^, respectively. These predictions may be tested with 
observations over a broad energy band. 

For the September 20 flare the total energy injection in the decay phase is ~ 6 x lO^^erg 
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( assuming a duration of ~ 10 mins), which is comparable to that in the impulsive phase 
(see § 2.3). The total energy injection in the decay phase of the August 12 flare (~ 10^® 
erg s~^ for a period of ~ 3 mins), on the other hand, is more than one order of magnitude 
lower than that injected during the impulsive phase, and the flare decay is dominated by the 
radiative cooling process although a suppression of conduction is still required. The Neupert 
effect therefore gives a reasonable description of this flare. These results, specially those for 
the September 20 flare, suggest that flares with a longer decay time require more energy 
injection in the decay phase and are not in line with the Neupert relation. 

Because it is believed that the energy of solar flares comes from the magnetic field, the 
low values of ^B^ inferred from the modeling (~ 130 and ~ 10 for the September 
20 and August 12 flares, respectively) suggest a ^ much less than 1 for the wave-particle 
interaction. This means much more efficient particle scattering than plasma heating by the 
waves. Once the wave modes responsible for these processes are identified, the observations 
can be used to constrain the wave energy density and spectrum (Schlickeiser & Miller 1998). 



4. Summary and Conclusions 

We have carried out imaging and spectroscopy studies of the evolution of 6 simple limb 
fiares, some of which reveal several interesting features of the LT source as enumerated below: 

• A distinct coronal LT source (Figures 1-3) is observed for each fiare. The source appears 
upon the onset of the fiare and exists until the fiare dies out. It is relatively stable in 
both shape and size and, in some cases, moves systematically toward high latitudes in 
the decay phase. 

• For a few fiares with their FPs completely occulted by the solar limb, the spectrum of 
the LT source can be studied without invoking imaging spectroscopy. Above 10 keV, 
the spectrum is dominated by a steep power-law during the preheating and impulsive 
phases (Fig. 6). A thermal component appears after the impulsive HXR peak and 
becomes more and more prominent as the fiare proceeds (Fig. 10). 

• The LT source spectrum can be fitted with a thermal model in the late decay phase. 
The cooling rate of the LT source is generally higher than the radiative cooling rate 
but much lower than the conductive cooling rate predicted by Spitzer or the so-called 
non-local conductivity, requiring a suppressed conductivity and/or continuous heating. 

• Imaging of this LT source indicates that the thermal hot plasma is confined within a 
small region near the LT and does not extend to the FP regions. 
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We have considered several possibilities to explain these observations and carried out a 
detailed investigation of the energy balance in the post impulsive loop: 

• A small filling factor can not explain the discrepancy between the high decay rate ex- 
pected theoretically and the observed cooling rate, when the effects of both conduction 
and radiation are included. It also can not produce the confinement of the LT source. 

• For loop models with a uniform conductivity because the conductive heat flux and 
the radiative cooling rate are very sensitive to the temperature, the loop must have a 
nearly uniform temperature profile to balance efficient radiative cooling near the cold 
FPs. This results in LT sources extending much closer to the FPs than those observed 
by RHESSI. 

• The suppression of conduction therefore must be confined to the LT region. Moreover 
for some fiares (i.e. the September 20 fiare) continuous energy input is also required 
to balance the rapid cooling in the loop legs. For fiares with relative longer decay time 
the energy input in the decay phase can exceed that in the impulsive phase, and most 
of these energy are radiated away near the FPs through the optical and UV emission 
fines. This prediction can be tested with broad band observations. The fluctuations 
in the energy decay rate and occasionally observed energy increase of some flares also 
suggest a continuous heating process. 

• We suggest that plasma waves or turbulence may account for both these aspects. The 
scattering of high energy electrons by waves can reduce the conduction. The particles 
can also gain energy from the waves, which would account for the heating. 

• Assuming the scattering time scale independent of the particle velocity and a resonant 
wave-particle coupling, which gives a heating rate inverse proportional to the square of 
the particle velocity, we showed that observations can be used to determine the size of 
the turbulence region and the degree of conduction suppression and plasma heating. 

This is a new result indicating that plasma waves or turbulence, which appears to play 
important roles in the acceleration of particles during the impulsive phase (Miller et al. 1987; 
Petrosian & Liu 2004), be generated continuously even during the decay phase, presumably 
at a lower level. We note that according to Petrosian & Liu (2004) a low level of turbulence 
means a slow plasma heating with little particles accelerated to high energies. The high gas 
density and temperature in the gradual phase also make the plasma heating more dominant. 
More detailed studies of wave-particle interactions, i.e. the scattering and acceleration of 
particles by PWT, will help us better understand the relevant phenomena, uncovering the 
energization mechanism of solar flares eventually. 
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Table 1: The ratios of tlie energy decay rates expected from the radiative, Spitzer conductive, 
and non-local conductive cooling to the observed one. The radiative ones are a bit smaller 
than unity while the conductive ones are much larger than one. The minimum energy 
decay rates due to the combined cooling of radiation and Spitzer conduction are obtained 
by adjusting the volume filling factor rj. The corresponding critical suppression factor 77^ 
is given in the last column. The fact that this minimum values are always larger than one 
shows that variations in rj can not explain the observations. 
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Fig. 1.— RHESSI PIXON images of the LT source (6-12 keV) and FP sources (at 51-57 keV 
and 34-39 keV, respectively) of the two hmb flares on September 20 and August 12. The thin 
contours (15% and 75% of the peak brightness of the image, which are adopted hereafter 
unless specified otherwise) and the gray scale are for the HXR peak (see light curves in figure 
10). The thick contours are for a period in the decay phase. The integration time for each 
image is ~ 20 seconds with the corresponding start times indicated in the figure. The dotted 
contours at the upper- right corner give the PSF at 6-12 keV. 
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Fig. 2. — Same as Figure 1 but for the partially occulted flare a and the corona event b on 
April 4. The grey scale is at 60-100 keV for the former, and for the latter the images are 
at 6-12 kcV with the grey scale for the HXR peak and the thin and thick contours for two 
periods in the decay phase as indicated. 
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Fig. 3. — Same as Figure 2 but for the two corona events on March 28 and April 30. 
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Fig. 4. — Left: Comparison of the simulated images with the observed image of the September 
20 2002 flare obtained with the PIXON algorithm. The gray scale shows the FPs (51-57 keV) 
at the HXR peak, and the red contours (15% and 75% at 6-12 keV) are for the LT with the 
time interval indicated. The thin black lines indicate the loop structure (curved) and the 
direction in perpendicular to it (straight). The simulated image contours are for the uniform 
loop model (blue contours), which is longer than the observed images, and for the model 
with a distinct LT region (green contours), which agree with the observations. Right: Same 
as the Left but for the August 12 2002 flare. 
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Fig. 5. — Left: Comparison of the observed and simulated image brightness profiles along 
(top) and in perpendicular to (bottom) the flare loop. The solid, dotted, dashed, and dot- 
dashed lines are for the observation, the non-uniform, uniform loop models, and the PSF, 
respectively. The LT is partially resolved in both directions. Right: Same as the Left but for 
the August 12 2002 flare. The LT is resolved along the loop, but the profile in perpendicular 
to the loop is consistent with the PSF. 
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Fig. 6. — The RHESSI total spectrum of the hmb flare of September 20 2002. The left 
panel is for the rising phase, fitted with a single power-law (lower), and for the HXR peak 
(upper), fitted with a thermal plus a power-law model, whose residual in units of the standard 
deviation is shown in the lower panel. The right panel is for the early decay phase (upper), 
fitted with a thermal plus a power-law model with the residuals in the lower panel, and for 
the late decay phase (lower), fitted with a pure thermal model. The corresponding time 
intervals are indicated in the figure and in the light curves by the arrows in the left panel of 
Figure 10. 
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Fig. 7. — Same as Figure 6 but for the August 12 flare. Here the peak time spectrum is 
also fitted with a single power-law model. The corresponding light curves are shown in right 
panel of Figure 10. 
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Fig. 8. — Same as Figure 7 but for the corona event b on April 4 2002. The corresponding 
hght curves are shown in right panel of Figure 11. 
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Fig. 9. — Spectra of the LT and FP sources of the September 20 (left panels) and August 12 
(right panels) flares during the HXR peaks. The solid hnes give the spectral fits to the FPs: 
a thermal plus a power-law model for the former and a power-law model for the latter. The 
LT spectrum is fitted by a thermal model as indicated by the dotted fine. For the August 12 
fiare the LT source does not yield a reliable spectrum in the range of figure and therefore is 
not shown. The ratio of the FP and LT fiuxes at lower energies ( < 20 keV) of the September 
20 flare are close to the dynamical range of the PIXON algorithm. The FP fluxes are thus 
not trustworthy. The lower panels show the relative error of the spectral fittings to the FPs. 
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Fig. 10. — RHESSI light curve and the evolution of model parameters of the hmb flares 
of September 20 2002 (left panel) and August 12 2002 (right panel). The top panels give 
the 6-12, 25-50, and 50-100 keV light curves with the higher energy ones shifted by the 
factors indicated in the figures for clarity. The arrows on the top indicate the times for the 
spectral fits shown in Figures 6-7. The second panels show the evolution of the emission 
measure EM (dashed line) and temperature T (solid line) for the thermal component, and 
the corresponding gas density n (dotted line) read from the observed source size. Here all 
quantities are given in cgs units and are scaled by the order of magnitude indicated in the 
subscript. The spectral indexes for either a single power- law fit or a broken power-law fit 
are also indicated with the dash-doted lines (gray for higher energy component and black 
for lower one for the broken power-law model). The third panel gives the observed energy 
decay rate (pluses) and those due to the radiative cooling (dot-dashed), Spitzer conduction 
(thick solid) and non-local conduction (thin solid). 
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Fig. 11. — Same as Figure 10 but for the partially occulted flare a (Left) and the occulted 
flare b (Right) on April 4 2002. 
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Fig. 12. — Same as Figure 10 but for the flare on March 28 2002 (Left) and the flare on April 
30 2002 (Right). The star signs indicate the energy increase rate during the periods. 
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Fig. 13. — Left: The temperature profiles (top panel) for several isobaric static loop models. 
The top dashed line is for a model without suppression of conduction and with a heat flux 
injected at the LT Fi^. For the rest profiles turbulence is assumed to distribute at the LT, and 
the ratio of scattering rates of the electrons by the waves and background particles t~^/t^^^^ 
is proportional to T^/^ and has a Gaussian form with a width w — 4'/5 and a peak value of 
So along the loop. The top panel shows how the temperature profile changes with 5*0, which 
is, respectively, 15, 9, 5, 2 for the solid lines from bottom to top. The model shown by the 
thick solid line is for the maximum possible conduction suppression (see text) and is then 
used for imaging simulation in Figure 4 and 5. The lower panel shows how the energetics 
changes with Sq. Here the total heating rate by the waves F^^ = J SUekBT < t~^^ > dV, the 
radiative cooling rate Sj-ad = J ^raddV^, the observed energy decay rate S and the conduction 
heat fiux at FPs Fout are designated by the thick sohd, dotted, dashed and thin sofid lines, 
(and by the plus, cross, diamond, and triangle signs for the model without a conduction 
suppression,) respectively. Right: Same as the Left hut for the fiare on August 12 2002. Here 
w — 6'/5, and in the upper panel -S'o = 4, 10, 18, 28 for the sohd lines from top to bottom. 
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Fig. 14. — Left: The energetics along the loop of the model with a maximum conduction 
suppression for the September 20 flare. The conductive heating flux is indicated by the thick 
solid line. The thin solid, dotted, and dashed lines give the heating, radiative cooling, and 
energy decay rates integrated from the LT, respectively. Note that more than half of the 
radiative energy is carried away near the FPs. Right: Same as the Left but for the flare 
on August 12 2002. Heating is less important here, and the observed energy decay can be 
attributed to cooling near the FPs. 



